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ABSTRACT 

The Kamioka Liquid-scintillator Anti-Neutrino Detector (KamLAND) 
has detected for the first time the disappearance of electron antineu- 
trinos from a terrestrial source at the 99.95% C.L.^ Interpreted in 
terms of neutrino oscillations i?e ^ i^x, the best fit to the KamLAND 
data gives a mixing angle sin^ 29 = 1.0 and a mass-squared difference 
Am^ = 6.9 X 10~^ eV^, in excellent agreement with the Large Mix- 
ing Angle (LMA) solution to the "solar neutrino problem".^ Assuming 
CPT invariance, this result excludes other solutions to the solar neu- 
trino problem at > 99.95% C.L. 



1 Reactor Antineutrino Experiments and Neu- 
trino Oscillation 



Nuclear reactors emit a calculable flux of electron antineutrinos (z^es) in all direc- 
tions. For standard particle propagation, one expects a detector located a distance 
L from the reactor to measure a flux that decreases as But if UgS are mas- 

sive, they may "oscillate" into undetectable flavors on the way to the detector, 
leading to an apparent "disappearance" of the z/gS. 

Neutrinos are produced and detected in weak interactions, which couple to the 
weak eigenstates ui, where I — e, /i, t. For massive neutrinos, the weak eigenstates 
may be expressed as a linear combination of three mass eigenstates Ui, i — 1,2,3, 
with mass mj: 

i^i^^UiiUi. (1) 

i 

is a 3 X 3 unitary "mixing" matrix and is analogous to the CKM matrix in the 
quark sector. As a neutrino propagates through vacuum with momentum p^, the 
phase of each mass eigenstate will change at different rates according to 

ui{t)^J2^ue-'^^\, (2) 

i 

2 

where Ei ^ py + At times t > the neutrino will be in a superposition of the 
weak eigenstates. The probability of detecting flavor t/;/ at a distance L k, t from 
the source is then 

2 2 

P{vi ^ vr) ^ReJ2 UiiUlUijUj.e'^^^ (3) 

In neutrino oscillation experiments, it is common to simplify to two neutrino 
flavors, say ly^ and u^, in which case Uu is parameterized by a single mixing angle 
9. In this picture, the probability of a neutrino emitted as to be detected as 
is then 

P {ve ^ i^e) = 1 - sin^ 29 sin^ ( ^ ) (4) 



osc 



where Lose — Here, Am^ = \ml — ml\ is taken to be the mass-squared 

difference between two "effective" mass eigenstates relevant to i/g <-> i/^ mixing. 
Note that P {lyg ^e) is a sinusoidal function of L, hence the name "neutrino 
oscillation" . 



The situation becomes slightly more complicated for neutrinos propagating 
through matter. As first recognized by Wolfenstein,'^ Mikheyev, and Smirnov,"^ 
while all three weak states participate in neutral-current interactions with normal 
matter, only electron neutrinos have additional charge-current interactions with 
electrons in the material being traversed. This results in an effective index-of- 
refraction for Vf, different from that for z/^ t- by 

An=H^/(0), (5) 

where A^e is the number-density of electrons in the material, and /(O) is the forward 
scattering amplitude. iSn changes the phase of the fe-component relative to the 
other components by 27r over a distance Lq given by 

where Gp\s, the Fermi constant. This phenomenon, called the MSW effect, alters 
the two flavor oscillation probability as follows: 

P {v, ^ i/e) = 1 - sin^ (2^^) sin^ (7) 
tan = tan 29(1 + ^ sec 29\ (8) 
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(9) 



For neutrinos created in and propagating out of the sun, Lq ~ 200 km <S Rsun, 
and matter effects are significant. The result is that Ve disappearance experi- 
ments that use the sun as a source, such as Homestake,^ GALLEX,^ SAGE,^ 
Kamiokande,^ Super-Kamiokande,^ and SNO,^° individually allow values of tan^ 9 
and Am^ over many orders of magnitude. The overlap of allowed regions from 
multiple experiments with different thresholds leaves a few small patches, with 
values around tan^^ Ri 1 and Am^ 10~^ eV^ highly favored. This region is 
called the "Large Mixing Angle" (LMA) MSW solution. Prior to KamLAND, 
another patch at lower Am^ ^ 10~^ eV^, known as the "LOW" MSW solution, 
survived at the 99.73% C.L.^° 

For neutrino experiments with artificial sources on the earth, matter effects 
are typically much less significant. This is because for matter with density similar 



to rock Lq ~ 10'^ km, larger than the radius of the earth. Thus for these sources 
the modifications to the oscillation parameters can be viewed as perturbations on 
the vacuum parameters, and for short enough oscillation lengths can be neglected. 

Artificial v sources for oscillation experiments typically take one of two forms: 
neutrino beams or nuclear reactors. These two sources comprise a highly comple- 
mentary experimental program in neutrino oscillation. Neutrino beams in gen- 
eral require very long baselines for good Am^ sensitivity, but the collimation of 
these sources lessens the severity of this obstacle considerably. Moreover, beams 
produce higher energy neutrinos and multiple fiavors, allowing for not only disap- 
pearance but also appearance measurements. Reactors, on the other hand, are a 
non-collimated source, and produce exclusively electron (anti) neutrinos, restrict- 
ing the experimenter to disappearance-only measurements and thus giving limited 
sensitivity to sin^ 1Q. However, due to the lower energies of the neutrinos (1-10 
MeV), reactor experiments are sensitive to very small Am^. 

Note that both neutrino beams and reactors are "laboratory-style" experiment, 
in which both the source and the detector are controlled. The physics of neutrino 
propagation can then be cleanly tested separately from the mechanisms involved 
in the production of the neutrinos. 

2 Calculating Reactor Antineutrino Fluxes 

More than 99.9% of z/gS emitted by nuclear reactors are produced in the decay 
chains of the fission products of only four isotopes: ^^^U, ^^^U, ^^^Pu, and ^^^Pu. 
The number of fissions of each of these isotopes over the data-taking period is 
estimated from a Monte Carlo simulation of the reactor core that tracks fuel 
burn- up and U/Pu production over reactor fuel cycles. The performance of such 
simulations has been verified by comparing the simulated fuel composition at the 
end of a fuel cycle with measurements of isotopic abundances in the actual spent 
fuel rods. 

The inputs to the simulation include the initial fuel composition and periodic 
measurements of the secondary calorimetric power, the pressure and flow rate in 
the primary cooling system, and various other operational parameters. However, 
the results of the simulations depend very weakly on most of the inputs. Accuracy 
is required only for the fuel composition and the power, on which the fission rates 
depend linearly. The power is measured rather precisely: reactors are regulated 



to operate at about a standard deviation or so below their rated power outputs. 
The economic incentive to produce as much power as possible pushes this error 
to be small; uncertainties < 1% are common. 

The emitted spectrum is then calculated by summing the spectrum emitted 
by each isotope weighted by the number of fissions of that isotope during the data 
taking period. For ^^^U, ^^^Pu, and ^"^^Pu, spectra are derived from /3-spectra 
measurements. ^^'^^ Extracting the z/g spectra from the /3-spectra is non-trivial; 
it involves complex bookkeeping of the decay chains for each daughter nucleus 
and summing the contributions. The spectra are normalized to the number of 
z/gS emitted per fission, and have uncertainties of a few percent. For ^^^U, no 
measurements are available, so we must rely on a calculation^^ of the Ve spectrum. 
The uncertainty in the calculation is ~10%, but since fissions from 238u j^iake up 
only ~10% of the signal, this results in a 1% uncertainty in the full z/g fiux. 

Reactor experiments typically detect z/gS via inverse-/3-decay, z/g + p — ^ n + e"'", 
because the 2-particle final state distinguishes the z/g from most other particles 
and dramatically reduces backgrounds. The positron carries away almost all of 
the energy of the incident neutrino, so that the detected positron spectrum is to 
a good approximation the incident neutrino spectrum folded with the detection 
cross section and shifted by a constant energy. The cross section has a threshold 
of 1.8 MeV, and has been calculated^"^ with uncertainties on the level of a percent. 

The most recent generation of reactor experiments, Palo Verde^^ and Chooz,^^ 
showed that uncertainties on the level of a few percent can be achieved for these 
calculations. The fact that z/g disappearance has not been detected in reactor 
experiments with baselines up to ~1 km^^^^° allows one to view these experiments 
as a verification the reactor fiux calculations to within a few percent over most of 
reactor neutrino spectrum. 

3 The KamLAND Experiment 

KamLAND uses the entire Japanese nuclear power industry as a 180 GWth long- 
baseline source. The experiment is located underneath Mt. Ikenoyama in Gifu 
prefecture in central Japan. 80% of the z/g fiux comes from reactors at baselines 
of 140-210 km. This range of baselines makes KamLAND particularly sensitive 
to values of Am^ corresponding to the LMA MSW solution to the solar neutrino 
problem. Matter effects in rock can be neglected at this distance scale. 



The results presented here are for 145.1 Uve-days of data taking between March 
and September 2002. Calculating the incident flux at KamLAND during this 
period requires summing the flux from every reactor in Japan. Reactor data, 
including instantaneous power and fuel burn-up, is provided by the Japanese 
power companies. The flux from reactors in Korea and the rest of the world is 
estimated from reported power generation, and makes up only a few percent of 
the overall flux. 

As stated earlier, KamLAND detects z/gS via inverse-/5-decay, Ue + P ^ n + e"*", 
with the energy of the positron Ee+ related to the energy of the incident by 
essentially a constant offset of 1.8 MeV. The positron quickly deposits this energy 
in KamLAND's liquid scintillator before annihilating with an electron, giving a 
prompt energy deposit of E^, — 0.8 MeV. Meanwhile, the neutron thermalizes and 
eventually captures on a proton, releasing a 2.2 MeV gamma in the process. The 
characteristic capture time of ~ 210 fis between the prompt and delayed events 
provides a powerful background reduction. 

KamLAND's liquid scintillator is composed of 80% dodecane, 20% pseudoc- 
umene, and 1.52 g/liter of PPO. One kton of scintillator is held inside a 6.5 m 
radius transparent balloon made of a 135 fim thick sandwich of nylon and ethylene 
vinyl alcohol copolymer fllms. 1879 17- and 20-inch photomultiplier tubes (PMTs) 
mounted on the inner surface of an 18 m stainless steel sphere view events inside 
the scintillator. For the results reported here, only the 1325 17-inch tubes are 
used, corresponding to 22% coverage. A buffer of mineral oil between the PMTs 
and the balloon provide buoyancy for the balloon and shielding from radioactivity 
in the PMT glass. Outside the steel sphere is a water cherenkov outer detector 
(OD) serving as a muon veto. A chimney and deck at the top of the detector pro- 
vide access for calibration devices to be deployed along the z-axis. A schematic 
of the detector is given in Figure 1. 

An event is recorded whenever 200 or more PMTs report a signal within 125 
ns. A delayed-event window is then opened for 1 ms with a lower threshold of 
120 tubes. The pulses are digitized by Analog Transient Waveform Digitizers 
(ATWDs). There are two ATWDs per PMT to reduce dead time during digitiza- 
tion. 

The visible energy and position of each event is reconstructed from the col- 
lected charge and the hit pattern. The reconstruction algorithms are developed 
and tested using uniformly distributed spallation products following muons and 
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Figure 1: Schematic diagram of the KamLAND detector. 

using cahbration data along the z-axis with ^®Ge, ^°Co, and ^^Zn gamma sources, 
and an AmBe neutron/gamma source. The cahbration data are also used to esti- 
mate the conversion from visible energy to kinematic energy, taking into account 
particle-dependent non-linear effects due to quenching and cherenkov light pro- 
duction. The systematics on the energy scale is shown in the upper panel of 
Figure 2. 

Antineutrino events are selected according to the following criteria: time cor- 
relation 0.5 /IS < At < 660 fj,s; vertex correlation AR < 1.6 m; prompt event 
energy Ep > 2.6 MeV; delayed event energy 1.8 MeV < < 2.6 MeV; spherical 
fiducial volume i? < 5 m; cylindrical radius p > 1.2 m. The cut on Ep elimi- 
nates Pe backgrounds from geological U/Th sources. The cylindrical radius cut 
removes backgrounds from thermometers deployed along the z-axis to monitor the 
scintillator. The total efficiency of these cuts was determined to be 78.3 ± 1.6%, 
and the residual accidental backgrounds, estimated using an off-coincidence time 
correlation window 20 /is < At < 20 s, is negligible at < 10~^ events per day. 
The size of the fiducial volume is determined by counting the ratio of uniformly 
distributed spallation products that pass the fiducial volume cut. The distribution 
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Figure 2: (a) Fractional error in energy reconstruction relative to the known source 
energies. The dashed line shows the systematic uncertainty in the energy scale, 
(b) Distribution in R"' of 2.2 MeV capture gammas from neutrons following muons. 
The ratio of events occurring inside the fiducial volume cut gives the fractional 
size of the fiducial volume. 



for spallation neutron capture gammas is given in the bottom panel of Figure 2. 

The correlated backgrounds passing the u^, cuts are produced by cosmogenic 
spallation. An overburden of 2700 m.w.e. reduces the muon rate to ~ 0.3 Hz. 
Muons arc identified by their large energy deposit and outer detector activity, and 
the muon track is reconstructed from the PMT hit pattern. A 2 ms veto following 
muons eliminates backgrounds due to spallation neutrons and short lifetime spal- 
lation products. Long lifetime spallation products that mimic the z/g signal, such 
as ^He and ^Li, are removed by a 2 second veto in a 3 m radius around the muon 
track, plus a 2 second veto over the entire volume following muons with very high 
energy deposits (>3 GeV). The residual correlated backgrounds, estimated from 
spallation event studies, time-since-last-muon distributions, and, for fast neutrons 
following muons not tagged by the OD, Monte Carlo simulations, is 1 ± 1 event 
in 145.1 days. 

Figure 3 shows the distribution of the z/g candidates' prompt and delayed 
energies. The delayed 2.2 MeV gamma cleanly separates Ue from the accidental 
coincidences clustered at small energies. The one event with Ed ^ 5 MeV is 
consistent with neutron capture on carbon; such events are not considered in this 
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Figure 3: Prompt and delayed energies of z/g candidates. 
Table 1: Estimated systematic uncertainties (%). 



Total LS mass 


2.1 


Reactor power 


2.0 


Fiducial mass ratio 


4.1 


Fuel composition 


1.0 


Energy threshold 


2.1 


Time lag 


0.28 


Efficiency of cuts 


2.1 


u spectra^^"^^ 


2.5 


Live time 


0.07 


Cross section-^^ 


0.2 


Total systematic error 






6.4% 



analysis at a small cost in efficiency with negligible uncertainty. 

The systematic uncertainties are listed in Table 1. The largest contribution 
to the systematics comes from the estimate of the fiducial mass ratio. The total 
uncertainty is 6.4%. 

For no oscillations, in 145.1 live-days we expect 86.8 ± 5.6 events with 1 ± 1 
of these coming from backgrounds. The observed number of events is 54, giving 
a ratio = 0.611± 0.085 (stat) ± 0.041 (syst). This ratio is plotted vs. 

baseline for KamLAND and previous reactor experiments in Figure 4. The shaded 
region corresponds to the LMA MSW solution to the solar neutrino problem, 
with which the reactor experiments are in excellent agreement. "Standard" Ue 
propagation is excluded at > 99.95% C.L. 
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Figure 4: Nobs/Nno osc for reactor z/g experiments. 



The positron energy spectrum is histogrammed in Figure 5. An unbinned max- 
imum hkehhood fit was performed and included terms for changes in the overall 
rate and shape of the spectrum at different points in Am^ — sin^ 26 parameter 
space. The best fit solution is plotted as the shaded histogram, and corresponds 
to Am^ = 6.9 X 10^^ eV^ and sin^26' = 1.0. The regions in parameter space 
allowed at 95% CL. are drawn in Figure 6. Also shown is the exclusion region 
in Am^ — sin^ 26 parameter space based on the KamLAND rate alone. Assuming 
CPT invariance (i.e. that v and v masses are identical), all neutrino oscillation 
solutions to the solar neutrino problem are excluded except LMA. The KamLAND 
result divides LMA into two regions of higher and lower Am^. The sensitivity in 
sin^ 26 is rather poor. Values of Am^ between the inclusion regions in general 
predict large distortions in the energy spectrum that were not observed and are 
hence disfavored. Constant suppression of the non-oscillated spectrum is consis- 
tent with the data at 53% CL. The inclusion region extends to high Am? where 
spectral distortions can not be distinguished with current statistics and energy 
resolution. An upper limit can be placed on Am^ using the non-observation of Ve 
disappearance by Palo Verde and Chooz. 
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Figure 5: Prompt energy spectrum. 
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Figure 6: 95% inclusion region from KamLAND rate + shape. 
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Figure 7: Time variation of the z/g flux at KamLAND. 

4 Future of the KamLAND Reactor Measure- 
ment 

Figure 7 shows the Ve flux at KamLAND from March 2002 through May 2003. 
In early 2003, Japanese utilities powered down their reactors for inspections and 
maintenance, resulting in a drop in the non-oscillated flux at KamLAND by about 
a factor of two. Such a drastic time variation will provide a more precise estimate 
of backgrounds, and depending on the "true" values of Am^ and sin^ 25, may 
also provide a time-varying spectral shape that will further restrict the inclusion 
regions in parameter space. 

The expected sensitivity of KamLAND for 5 years of data taking is shown 
in Figure 8 for two sets of parameters expected from KamLAND's first results. 
KamLAND stands to provide exquisite sensitivity in Am?. However, further 
restrictions on sin^ 26 will probably have to come from future solar neutrino ex- 
perimental results. 
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Figure 8: Expected KamLAND sensitivity for 5 years of data. 

5 Future of Reactor Experiments 

With the latest generation of solar and reactor experiments dawns a new age in 
precision neutrino oscillation physics. The next step in this exciting field will be to 
improve on current measurements and to measure some of the remaining unknown 
parameters in full three-flavor mixing. In addition to the three masses m^, there 
are 6 free parameters in the mixing matrix Uu. It is possible to parameterize Uh 
as follows: 

U = UlmA X Uatmospheric X f^l3 X Upp. (10) 
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Uatmospheric and the value of Am^g has been measured by the Super-Kamiokande 
experiment. Ufjp and the overall mass scale will require kinematical and neutrino- 
less double-/? decay measurements. KamLAND is playing a central role in deter- 
mining Ulma as well as the value of Amf2- Reactor experiments stand to play 
another role in neutrino oscillation, this time in the measurement of [/13. 

Because Am^2 Am23, it must be the case that Am^g ^ Amgg. A new 
generation of reactor experiments has been proposed to search for z/g disappearance 
at baselines of 1 km corresponding to this value of Am^. To improve on the 
mixing angle sensitivity achieved by Palo Verde and Chooz. proposals for reactor 

experiments include a large detector to reduce the statistical error, and also a 
second detector positioned very close (~100 m) to the reactor. The near detector 
would precisely measure the incident flux, allowing many of the flux calculation 
systematics to drop out. This also requires that the detectors be made identical 
and/or moveable. Sensitivity down to sin^ 2^13 ^ 10^^ seems within grasp. Such 
experiments were first discussed by Mikaelyan and Sinev^^; for a comprehensive 
list of references on this topic, please see the web site compiled by Heeger.^^ 



6 Conclusion 

KamLAND has observed, for the first time, disappearance of electron antineu- 
trinos in a lab oratory- style experiment. Assuming CPT invariance, this result 
excludes solar neutrino oscillation solutions except LMA at > 99.95% C.L. Recall 
that the LMA solution to the solar neutrino problem is a region in oscillation 
parameter space allowed from the overlap of many experimental results using 
different techniques and with different thresholds. Moreover, LMA is an MSW 
solution for neutrinos, in which matter effects inside the sun drive the oscillations. 
It is significant that KamLAND, which uses yet another detection technique and 
is sensitive to vacuum oscillations of antmeutrinos, gives oscillation parameters 



in agreement with LMA. These experiments are different in so many aspects, yet 
the physics of neutrino oscillation ties them together in a beautifully consistent 
theoretical framework. 



7 Acknowledgements 

The KamLAND experiment is supported by the Center of Excellence program 
of the Japanese Ministry of Education, Culture, Sports, Science and Technol- 
ogy, and the United States Department of Energy. The reactor data are provided 
courtesy of the following electric associations in Japan: Hokkaido, Tohoku, Tokyo, 
Hokuriku, Chubu, Kansai, Chugoku, Shikoku, and Kyushu Electric Power Compa- 
nies, Japan Atomic Power Co., and Japan Nuclear Cycle Development Institute. 
Kamioka Mining and Smelting Company provided services for activities in the 
mine. 

References 

[1] K. Eguchi et ai, Phys. Rev. Lett. 90, 021802 (2003). 

[2] J.N. Bahcall and R. Davis, Science 191, 264 (1976); J.N. Bahcall, Neutrino 
Astrophysics (Cambridge University Press, Cambridge, UK, 1989); J.N. Bah- 
call, Astrophys. J 467, 475 (1996). 

[3] L. Wolfenstein, Phys. Rev. D 17, 2369 (1978). 

[4] S.P. Mikheev and A.Yu. Smirnov, Sov. J. Nucl. Phys. 42, 913 (1985). 

[5] B.T. Cleveland et al, Astrophys. J 496, 505 (1998). 

[6] GALLEX Collaboration, Phys. Lett. B 447, 127 (1999). 

[7] J.N. Abdurashitov et ai, Phys. Rev. C 60, 055801 (1999). 

[8] Y. Fukuda et al, Phys. Rev. Lett. 77, 1683 (1996). 

[9] Super-Kamiokande Collaboration, Phys. Lett. B 539, 179 (2002). 

[10] Q.R. Ahmad et al, Phys. Rev. Lett. 89, 011302 (2002). 

[11] A.A. Hahn et al, Phys. Lett. B 218, 365 (1989). 

[12] K. Schreckenbach et al, Phys. Lett. B 160, 325 (1985). 

[13] P. Vogel et ai, Phys. Rev. C 24, 1543 (1981). 



[14] P. Vogel and J. F. Beacom, Phys. Rev. D 60, 053003 (1999); radiative cor- 
rection from A. Kurylov, et al., hep-ph/0211306. 

[15] F. Boehm et al, Phys. Rev. D 64, 112001 (2001). 

[16] M. ApoUonio et al, EPJC 27, 331 (2003). 

[17] H. Kwon et al, Phys. Rev. D 24, 1097 (1981). 

[18] G. Zacek et al, Phys. Rev. D 34, 2621 (1986). 

[19] G.S. Vidyakin et al, JETP 59, 390 (1994). 

[20] B. Achkar et al, NucL Phys. B 434, 503 (1995); B. Achkar et al, Phys. Lett. 
B 374, 243 (1996). 

[21] Super-Kamiokande Collaboration, Phys. Lett. B 467, 185 (1999). 

[22] L.A. Mikaelyan and V.V. Sinev, Phys. Atom. Nucl. 63, 1002 (2001). 

[23] http://thetal3.lbl.gov/references.html 



